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Lattice Investigations of the Running Coupling 

p. Weisz'' 

^Max-Planck-Institut fiir Physik, Fohringer Ring 6, D-80805 Miinchen, Germany 

A review of investigations of running couplings using lattice techniques is given. This includes i) studies of the 
running of particular non-perturbatively defined renormalized couplings in pure gauge theories over a range of 
energies, and ii) how estimates of aj^{mz) in lattice QCD are presently obtained. 



1. INTRODUCTION 

A widely accepted property of QCD is asymp- 
totic freedom (AF), which allows the computa- 
tion of certain aspects of high energy E processes 
as a perturbative series in a coupling which de- 
pends on E and goes to zero as £" — cxd. There 
are infinitely many physically admissible non- 
perturbative definitions of running couplings. At 
low energies (LE) these can have completely dif- 
ferent behaviors from one another, but at high 
energy (HE) they behave similarly and can be re- 
lated to one another using perturbation theory 
(PT). 

The currently most widely used scheme for the 
analysis of experimental data in HE physics is 
the MS scheme [1] using dimensional regulariza- 
tion. However it is important to bear in mind 
that (Xj^ is only operationally defined in PT. Its 
extraction from experimental (or numerical) data 
is not entirely straightforward. Obviously it re- 
quires the dropping of higher order terms which 
have not been computed. But more awkward is 
the fact that the nature of the perturbative se- 
ries is unknown and there are principle uncertain- 
ties [2] due to the neglect of non-perturbative ef- 
fects; these are typically 0{A^j^/ E^) (where Aj^, 
the so-called lambda-parameter, is a characteris- 
tic hadronic scale), only a priori negligible at re- 
ally high energies. For example to take the uncer- 
tainty of such effects into account it has been sug- 
gested [3] that the error in aj^(m^) now quoted 
0.36(3) should be doubled (yielding 0.122(7) for 
aj^{mz) instead of 0.122(4)). 

Values of (Xj^ extracted from various exper- 
iments in an energy range up to the Z-peak 



are qualitatively consistent with the running pre- 
dicted by the renormalization group (RG) equa- 
tions using perturbative /3-functions, but the sta- 
tistical and systematic error bars are too big to 
make definite statements. A compilation of values 
for aj^{mz) from various experiments was made 
last year by Webber [4] . There is tendency for the 
central value extracted from low energy (and run 
to niz using RG) to be slightly lower than that 
obtained directly from the LEP collider at the 
Z-peak. An often posed question in present phe- 
nomenology [5] is whether this is due to underes- 
timating the systematic errors referred to above 
or whether the potential discrepancy should be 
taken seriously and interpreted as a manifesta- 
tion of new physics. 

A definition of QCD requires a non-perturba- 
tive regularization e.g. the lattice. The possibil- 
ity to compute running couplings in asymptoti- 
cally free theories through numerical simulations 
of the corresponding lattice theories has recently 
attracted much attention [6]- [33]. There are two 
main directions of research: in the first one ac- 
tually studies the running of a coupling over a 
range of energies; in the second one uses certain 
measurements (see sect. 6) to estimate (Xj^ at rel- 
atively low energies and then reaches a value of 
aj^{mz) assuming applicability of perturbative 
RG running. 

The definitions of the coupling which have been 
investigated on the lattice are given in sect. 2. If 
these are measurable to sufficient accuracy, one 
can actually test if the behavior is really that 
predicted by the RG with the perturbative /3- 
functions. There is no rigorous proof that this 
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must be the case, and some authors [34] favor 
the scenario that although the continuum limit 
of lattice QCD may exist and describe the LE 
physics properly, the HE behavior will not be that 
described by AF (i.e. typical couplings will not 
approach zero at HE). Such scenarios can practi- 
cally not be distinguished experimentally, but the 
question can be addressed on the lattice. 

Taking the conventional scenario as a working 
hypothesis, the determination of aj^{mz) with 
a reliable error estimate requires various basic in- 
gredients: 1) an accurate determination of a ref- 
erence energy scale (see sect. 3.); 2) the introduc- 
tion of a non-perturbatively defined renormalized 
coupling which can be measured accurately on 
the lattice over a wide range of energies; and 3) a 
careful treatment of the continuum limit. Hav- 
ing these one can see at which energy the ex- 
pected perturbative behavior actually sets in, and 
thereon one can confidently run to even higher 
energies e.g. mz, where one converts to the MS- 
scheme. In order to avoid introducing an addi- 
tional error > 1% at this final stage one needs 
the knowledge of the relation to (Xj^ to 2-loops 
(see sect. 4). 

The program above has so far only been achie- 
ved for pure gauge theories [7]-[13] (see sect. 5). 
To obtain preliminary estimates of aj^{mz) in 
QCD a shortcut has been proposed which involves 
the dropping of steps 2) and 3), and using instead 
a modified bare PT as explained in sect. 6, to es- 
timate (Xj^ at some lower energy and then using 
perturbative running to niz- It turns out in pure 
gauge theory that this procedure reproduces the 
results of the more detailed investigation, but one 
would not be able to say this a priori since the sys- 
tematic errors due to cutoff effects are not really 
controlled m the shortcut approach. 

In QCD the method above has been invoked 
[21]- [33] producing results (already quoted in the 
Particle Data Tables) consistent (at present) with 
experimental determinations. Such computations 
have been made in the quenched approximation 
and in full QCD with nj =2; these are reviewed 
in sect. 7. In sect. 8 we describe the status of the 
Alpha-collaboration project to study the actual 
running. Finally in sect. 9 we briefiy discuss the 
treatment of the heavy quark thresholds. 



2. DEFINITIONS OF COUPLINGS 

2.1. Couplings from the potential 

Various definitions of the coupling involve the 
static potential V{r) [14]-[19]. The first is 

amiq)=C^'^r^Fir), q=^/r, (1) 

where F[r) = AV{r)IAr is the force (up to a sign) 
and Cp = [N"^ — l)/2#. Alternatively one can 
consider the Fourier transform V{q) of the po- 
tential and define 

<Xv{q) = -^C^\'V{q). (2) 

Yet another recent suggestion [19] is to define cou- 
plings by using RxT Wilson loops with R/T = x 
fixed instead of taking the T ^ co limit. 

With these definitions it is difficult to get a{q) 
in a large momentum range, despite the fact that 
the quality of data on Wilson and Polyakov loops 
is now very good. 

2.2. From the 3-gluon vertex 

Couplings defined from the 3-point vertices 
were previously investigated in PT and consid- 
ered (despite gauge dependence) "more physical 
in some respects" than a-^. Recently a project 
has been initiated [20] to measure a renormalized 
coupling based on the triple gluon vertex G'-^-' (in 
the Landau gauge). Again the momentum range 
which can be covered is very small at present. 

2.3. Couplings running with the volume 

The first coupling which runs with the box size 
to be considered in gauge theories, is based on the 
Schrodinger functional (SF) [7] Z{L,ri). This is 
the partition function on a finite volume L x 
with periodic boundary conditions in the spatial 
directions and Dirichlet bc's in the 'time' direc- 
tion. The fields at the boundaries (xq = 0,L) 
have to be chosen appropriately (see [7]) and de- 
pend on a parameter rj. The coupling is then 
defined as a response of the system to changing 
the bc's 

d 

asF(g)"^ oc — lnZ(_L,?y)|^=o, q=l/L. (3) 
This coupling vanishes exponentially as _L — cxd. 
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A second definition is based on tlie ratio of cor- 
relations of Polyakov loops [11]: 



(4) 



where the gauge field satisfies periodic bc's in the 
z-direction but twisted bc's in the x, y directions. 
At low energies this coupling goes to a constant. 

These couplings seem at first rather exotic, but 
at HE they can be matched to infinite volume 
physical couplings using PT. The basic advan- 
tages of couplings which run with the box size L 
were explained in detail in ref.[9]. Firstly they 
can be accurately measured over a wide range 
of energies. Secondly, a very important aspect, 
is that there is a well-controlled extrapolation to 
the continuum limit. 

2.4. Generalities 

All the couplings above can be defined in the 
theories including fermions. Note, as mentioned 
previously, they all behave asymptotically the 
same at HE but very differently in the LE (non- 
perturbative) regime. Given one coupling one can 
trivially define from it another valid coupling hav- 
ing different LE behavior. An illuminating ex- 
ample was given by Michael last year [16]; for 
pure gauge theories consider the coupling aFil) = 
Cp^r^(F(r) — K), where K = F(oo) is the string 
tension. At HE aFil) and aqq(g) have the same 
perturbative expansion and differ only by a non- 
perturbative term cx exp[— l/47r&oaqq(9)] ^ which 
becomes negligible (< 1%) for q lOGeV. 

3. SETTING THE SCALE 

To set the energy scale one needs a physical 
infinite volume observable go (with mass dimen- 
sion). At a given bare coupling ag, a measure- 
ment of ago then gives the lattice spacing in units 
of go (or in units GeV""'^ if go is assigned its physi- 
cal value) . The dependence on the observable will 
not be considerable provided one is in the scaling 
region where cutoff effects are small. The working 
hypothesis, motivated by the work of Symanzik 
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[35], is that physical quantities approach their 
continuum limit with power corrections in a. 

The problem is to find a scale which can be 
measured accurately since errors Sqo in a scale go 
give errors Sa/a ~ 0.2Sqo/qo in « at niz i.e. one 
needs the scale to < 5% to avoid introducing an 
error > 1%. 

For the pure gauge theory one can consider for 
example a glueball mass or the string tension K . 
Starting from the potential V a better parameter 
is Sommer's scale rg [36]: 



rlF{ro) = 1.65. 



(5) 



It is easier to extract from the data than K 
and its definition extends to the theory with 
fermions. Actually rg is not obtained directly 
from experiment; the value 1.65 above is chosen 
since the phenomenological potential describing 
the spectrum of charmonium has this value for 
r-o ~ 0.5fm. For N = 3, ^Kro ~ 1.2 [36]. 

One could contemplate eventually using some 
spectral value nip or mp . But at present this is 
rather difficult because of the required quark mass 



extrapolation m„ 



0, and in the case of nip one 
needs to take into account resonance effects. Al- 
ternatively one could take meson decay constants 
e.g. /,r, but this involves careful non-perturbative 
evaluation of the axial current normalization con- 
stant. 

Setting the scale with mass splittings in heavy 
quarkonium [23,21] IP — 1S,2S — IS has vari- 
ous advantages. Firstly they empirically depend 
only weakly on the heavy quark mass and sec- 
ondly they are expected to be rather insensitive 
to very light quark effects. The main problem is 
that in presently conceivable simulations of rel- 
ativistic fermions rrica ~ 1, so that cutoff effects 
are potentially large. The way to overcome this is 
to use an appropriate effective lagrangian for the 
heavy quark system e.g. a non-relativistic treat- 
ment [28,21] (see also sect. 7). 

4. SCHEME DEPENDENCE & BOOSTS 

4.1. Renormalized coupling relationships 

As mentioned before, at HE different schemes 
are related by PT. For example through the series 

"ms(«9) = "x(g)+cf (s) ax(g)^+cf (s) ax(g)^+--(6) 
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Table 1 

Coefficients ci in eq.(6) for massless quarlss 



Scheme X (1) 



MS 




47r&o(ln 47r — 7^;) 


99 




(-31#+ 10n/)/367r 


V 




Cf (1) - 87TbojE 


SF, N = 


2,n/ = 


0.94327(4) 


SF, N = 


3 


1. 256+0. 023n/ 


TP, N = 


2,n/ = 


-0.5584(4) 



one can convert from the X-scheme to the MS 
scheme [1]. Since a ~ 0.1 at the Z-peak, it is 
necessary to know the coefficient if we want to 
obtain results for aj^{mz) to a precision of 1%. 

The s-dependence of the coefficients is fixed by 
the RG equation: 

cf(s) = -8^&o Ins + cf(l), (7) 

cf(s) = cf(s)2 -32^2^1 Ins + cf(l) - cf (1)2. (8) 

The ci have been computed for all the schemes 
mentioned in sect. 2 and are given in Table. 1. 

The reason for introducing the scale factor 
("boost") s in eq.(6) is the intuition that before 
making practical use of the series, s should be 
fixed so that the low-order terms are reasonably 
small (as far as possible). One can e.g. chose the 
value of s = So such that ci(s) vanishes; this is 
just given by the ratio of the A-parameters ^ 

cX(,X, = „,„„X^exp{|ili} = M. 

For the special case of the MS-scheme we have 
"ms(«o*^9) = "Ms(g) (10) 

to all orders. 

Apart from this rather trivial example, at 
present the only known 2-loop parameter C2 is 
for the scheme AT = SF for # = 2, ri/ = [37,38]. 

4.2. Expansions in the bare parameter 

To get the relationships between the different 
schemes it is often easier (albeit not necessary) to 

^Note that the value of C2(s^), if positive, is (at least for 
N = 2,3,nf = 0) very close to the minimum value of C2, 
since cf{s) > cf{s^) - (27r6i/6o)^ for all s. 



Table 2 

Values of perturbative coefficients 



Coefficient Value 



ki 


2.135730074078457(2) 


k2 


-2.8626215972(6) 


k3 


1.24911585(3) 


ai 


0.0063537 


0-2 


-0.0181240 


as 


0.0185223 



first expand each coupling in powers of the bare 
coupling ao = gl/4:Tr, 

ax(q) = ao + df(qa) + df(qa) aj^ + . . . (11) 

In the continuum limit the bare coupling and the 
lattice spacing a are simultaneously taken to zero 
in such a way that a is an exponentially vanishing 
function of ao- In this context it is consistent to 
drop all terms of order ,p > contributing to 
the coefficients d„, which then reduce to polyno- 
mials in In(ga) as determined by the RG equation 
cf eqs.(7,8). When we combine the relations for 
different schemes to eliminate ao the dependence 
on a also drops out. 

The coefficients d^[l) depend of course on the 
lattice action used. Here we consider only a 4- 
dimensional hypercubic lattice and for the pure 
gauge SU(#) theory the standard Wilson action. 
The coefficients d^ for all schemes can be ob- 
tained from the knowledge of cf and (if ^ which 
was computed nearly 15 years ago by A. and P. 
Hasenfratz [39]: 

rff(l) = rfr-cf(l), (12) 

dril) = -^+kiN, (13) 

with ki given in Table 2. 

The first 2-loop coefficient d2 to be obtained 
was for the SF-scheme for the case of SU(2), by 
Narayanan and Wolff [37] 

rff(l) -rf^P(l)2 = 2.5375(4). (14) 

An estimate for d^^[l) has been obtained by 
Ellis and Martinelli [40,41] on the basis of "tad- 
pole dominance" (as first proposed by Parisi [42]). 
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But a complete computation was only presented 
recently by Liischer and myself [38]. Due to the 
complications of lattice PT we decided to use the 
background field technique [43,44], which has pre- 
viously been applied to compute the 1-loop coef- 
ficient df^(s) [45]. 

In ref.[46] the background field method on the 
lattice is described: this discussion leans heavily 
on the work of Kluberg-Stern and Zuber [44] in 
the continuum, and on the proof of renormaliz- 
ability of lattice PT by Reisz [47] . The applica- 
tion to the computation of was given in [38]. 
For the evaluation of lattice Feynman integrals 
in the continuum limit, we employed a newly de- 
veloped position space technique [48]. A crucial 
ingredient of this is an observation of Vohwinkel 
that the free boson (lattice) propagator can be 
computed recursively algebraically in terms of its 
values at the corners of the unit cube at the ori- 
gin. 

To determine d^^[l) the background field pro- 
pagator must be worked out to 2-loops on the lat- 
tice, and matched with the corresponding quan- 
tity in the continuum theory using the MS scheme 
of dimensional regularization which has been pre- 
viously computed by Ellis [41] ^. 

The final result is 

drii)-dr{ir = ^+k,+k,N\ (i5) 

where ^2, ^3 are given in Table 2. 

4.3. The coefficients for the Plaquette 

The (normalized) plaquette expectation value 
P [49] is now known to 3-loops [50]: 

_P = 1 - piao -P2ao - Pstto + • • • , (16) 
with 



Pi = Cftt, (17) 
P'^ = -f^F# {^2k, + ^ + ^) , (18) 




^and checked by A. van de Ven and M. Liischer. 



with a„ given in Table 2, and ^2 is the same con- 
stant as in eq.(15). The coefficients a„ have been 
checked using MC methods [51]. These compu- 
tations (and those of the previous subsection) il- 
lustrate the rapid growth of the application of 
algebraic and numeric computer methods for the 
evaluation of lattice Feynman diagrams. 

5. INVESTIGATING THE RUNNING 

5.1. The SF- and TP-couplings 

Last year Guagnelli [13] presented the results 
of the Alpha-collaboration study of the universal- 
ity of the continuum limit in pure SU(2) Yang- 
Mills theory. The two couplings asF and axp 
(see eqs.(3,4)) were accurately determined [12] 
over a large range of momenta given in units 
of Sommer's scale rg [36]. A crucial feature to 
appreciate in these analyses is that the results 
are extrapolated to the continuum limit, using a 
multi-lattice approach. The method has been dis- 
cussed in detail in many places (see e.g. [10]). The 
idea is that a sequence of lattices with parameters 
(j3 = PijL/a = Li/a) are determined such that 
say the coupling as-p{L) remains fixed, and then 
for the same parameters as-p{2L) is measured and 
the resulting sequence is extrapolated to the con- 
tinuum limit. Here, (as referred to in sect. 3), the 
working hypothesis is that cutoff effects fall as 
powers in a/L. It was found that the data was 
certainly consistent with such a behavior. The 
procedure is then repeated for as-p{L) fixed at 
the previously extrapolated value of as-p{2L), in 
order to obtain the running. 

For a sequence of lattices with fixed value of 
0!sf{L) we also measured the corresponding val- 
ues of aTp(-^)- These approached a definite value 
(again in the expected way) providing additional 
evidence for the existence of the continuum limit. 

A remarkable property of both couplings (in 
the continuum limit) is that their evolution in the 
momentum range covered by the numerical simu- 
lations is well described by the 2-loop approxima- 
tion to the Callan-Symanzik /3-function [8,12] 

perfect fit of the SF-data was in fact obtained by 
adding an effective 3-loop term with coefficient h^^ = 

0.35(f2)/(47r)3. With cf^ and the Icnown value of 6^^ 
[52] we may now calculate the exact coefficient = 



6 



only axp at the lowest energy shows a slight de- 
viation. We also found that the perturbative re- 
lation between the couplings (unfortunately only 
known presently to 1-loop) works better if one 
chooses the corresponding boost s (cf sect. 4) to 
be equal to the appropriate ratio of A-parameters 
rather than just setting s = 1. These features can 
be seen in Fig.l. 




Figure 1. A comparison of (/^p(£') (circles) and 
(/gp(AsF-E'/ATp) (boxes). The dotted line is 2- 
loop evolution. 

Since asF runs with the expected RG behavior 
in the HE range measured, it is reasonable to con- 
tinue the integration up to even higher energies 
and there convert to the MS-scheme. If we choose 
s = Aj^/AsF, (which proved successful above), 
then at q = 200/ro (which corresponds to about 
80 GeV in physical units), the result using the 
formulae only to 1-loop is [12] (SU(2)) 



aTjs(200/ro) = 0.1288(15)(21), 



(20) 



where first error is the total statistical error, while 
the second is ~ a|p (indicating the expected or- 
der of magnitude of the 2-loop correction). With 
the 2-loop coefficients eqs.(14,15) at hand we can 
now do better. In fact one finds a very small 2- 
loop coefficient C2^(so^) = 0.0576(5), which has 
practically no infiuence on the central value of 



There is no reason to expect that 
is equal to , but it is reassuring that both coefficients 



0.1797(3)/(47r) 



Table 3 

Evaluations of a- 
Method 



j^s(37/ro) for SU(3), ri/ 



SF,l-lp,s= 1 

SF,l-lp,S: 

aqq,/3 = 6.5;2-lp ev. 

eq.(28),/3 = 6.5;l-lp 
eq.(28),/3 = 6.5;2-lp 
eqs.(29-31),/3= 6.0 



0.1108(23)(14) 
0.1141(24)(14) 

0.1180(21) 

0.1121(14) 
0.1134(2) 
0.1113(11) 



+0 ■ 

-60 , 



eqs.(29-31),/3= 6.4,6.5 0.1126(10) 



[q) and the result eq.(20) is modified to 



aT7s(200/ro) = 0.1289(15)(3). 



(21) 



The main impact of the 2-loop correction is that 
the estimated error from the neglected higher or- 
der terms (now ~ alp) has been reduced from a 
few percent to a fraction of a percent and is now 
less than the statistical error. 

The SF-coupling has also been measured in the 
SU(3) theory [10]. Again the expected running 
was observed and the results are qualitatively 
similar to the case of SU(2). Unfortunately the 



3 is not yet available 



result for for N 
formulae are presently restricted to 1-loop. The 
result for (Xj^ from the data point at highest en- 
ergy is given in Table 3 for two choices of the 
boost s. 

5.2. Couplings from the potential 

The data on Wilson loop expectation values 
and on their correlations now available is of very 
high quality, both in SU(2) and SU(3) Yang-Mills 
theory, thanks largely to the Wuppertal group 
[17] and UKQCD [14,15], so that one now has 
information on the force and on the shape of the 
fiux tube [53], over a range of distances. If the 
force is defined as in eqs. (2.5-2.7) of ref.[36] then 
the cutoff effects seem to be rather small. The 
behavior of the resulting running coupling aqq at 
short distances is consistent with the expected 
perturbative RG behavior. Instead of showing 
ftqq, I would here like to show Fig. 2, the force it- 
self as a function of distance in units of ro, which 
emphasizes the non-perturbative aspects. Fig. 2 



have the same sign and order of magnitude. 



the computation is under way by A. Bode 
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was produced by Liischer from data (in Table 3 
of ref.[36]) for /3- values 2.5-2.7. The dotted line is 
the curve (r~^ — r[^^)7r/12 + 1.65r[^^. One indeed 
observes remarkable evidence for the presence of 
Liischer's universal 7r/(12r^) force [54]! Also it 
is rather surprising that this formula, which is 
an asymptotic formula for the non-perturbative 
large r regime, actually describes the data quite 
well up to quite HE - yet another example of pre- 
cociousness (but this time from the other end). 



of the ansatz, and within such a procedure it is 
hard to make reliable error estimates. 

5.3. The MOM-coupling 

A preliminary study of the running of the 
MOM coupling in SU(3) Yang-Mills theory has 
recently appeared [20]. Unfortunately the mo- 
mentum range covered at present is rather lim- 
ited, but the authors claim their results are in 
satisfactory agreement with those obtained by the 
other methods. 



6. USES OF BARE PT 



O IV) CO 



very 



Figure 2. Force in SU(2) pure gauge Iheory, 

The force in the SU(3) [55] thgory 4ook 
much the same as in Fig. 2; the co^e§]3 indint plot 
(updating that in ref.[55]) can be obi aine(| from 
Wittig. As a typical result for aj^,ji^aEe th| data 
for aqq[q) at /3 = 6.5, which covers tl-e rai|ge up 
to g ~ 7/ro. Then using the 1-loop perturiative 
relation (with s = 1) to go to the ItfS-s^heme 
and 2-loop evolution from thereon one obtaiSis the 
third entry in Table 3, in satisfacto_^ -agrefement 
with the SF-value. Ka 

There are many papers which rely on global 
fitting procedures to extract the coupling. For 
example Klassen [18] proposes to obtain ay by 
fitting the lattice data for the potential with a 
simple ansatz which assumes a particular form for 
the cutoff effects. The the expected continuum 
behavior in IR and UV is also built in and the 
associated parameters are obtained by optimizing 
the fits. However attractive a particular ansatz 
may be, the results depend sensitively on the form 



At first it seems that the series eq.(ll) is of lit- 
tle practical use since it only applies if ag In(ga) 
is small. But, as explained in [12], a similar sit- 
uation is encountered in the continuum theory, 
when ax{g) is expanded in powers of ax{g') with 
q' ^ q. If we now choose q' proportional to 1/a 
*he series may be combined to give 

I I I I I I I I I I I I I I I I I I I I I I I I I I 11 Tl 

cgtd,s.h)%ao + d^{s)al + d^{s)al + .... (22) 

. a iiJ 

— A'*^ < 

It shotlcfebe emphasized that in eq.(22) ax(s/a) 
is defir etf in the continuum theory. In particu- 
lar the momentum s/ a must be given in units of 
some physical reference scale (see sect. 3) for the 
relatio i to become useful. 

The formula above forms the basis of a more 
direct way to obtain the value of the running cou- 
pling a; say the Z-peak pioneered by the group at 
FNAL [21,56]. The idea is 1) to measure a phys- 
ical sccle at a given bare coupling to get a value 
of a[a[ ); 2) use a formula of the type of eq.(22) 



to get i value for ftf^ at q = s/a and finally 3) 
run to niz using the perturbative evolution. 

The question is, what is the best way to pro- 
ceed explicitly? If we choose for N = 3 e.g. 
s = Aj^/Aiat ~ 28.8 (as suggested by previous 
experience), one gets d^^[s) ~ 4.45. In view of 
the size of this 2-loop correction and the large 
scale factor one may hesitate to make practical 
use of this form of the expansion. The problem 
with the scale factor disappears if instead of ag 
one uses the "improved" bare coupling [42,56] 



ao = ao/P 



(23) 
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as the expansion parameter ®. Choosing again 
the scale factor so that the 1-loop term vanishes 
we end up with 

aj^(2.63285/a) = 5o + 2.18505(1) x + . . . .(24) 

This expansion looks more comfortable although 
in the relevant range of bare couplings, (/q — 1, 
the 2-loop correction is still 3-4%. 

Lepage and Mackenzie [56] now propose using 

an = -Pi^lnP. (25) 

Since the expansion of _P in ag is known to 3-loops 
the corresponding expansion of «□ is known to 2- 
loops. e.g. for N = 3,nf = we have 

«□ =ao + 3.373ag + 17.70aj^+... (26) 

At /3 = 6.5 this accounts for only 92% of the mea- 
sured value, of which the last contribution gives 
7%, thus this series seems poorly convergent. 

However it is observed [56] that physical cou- 
plings have better perturbative series in parame- 
ters such as an , which is of course also not a phys- 
ical coupling but is expressed purely in terms of a 
measured quantity and as such may include some 
'generic' non-perturbative effects. Going back for 
the moment to the case N = 2 one obtains 

asF(a"^) = an+ 0.731a^ + 1.125a^ + .... (27) 

If one now uses the UKQCD value «□ = 0.1479(3) 
at 13 = 2.85, one gets asF(20.6/ro) = 0.1677(6), 
which compares favorably to the value 0.1693(23) 
which is obtained from the extrapolation of the 
Alpha-collaboration data to that point. System- 
atic deviations show up however when the proce- 
dure is repeated at smaller /3. 
Returning to # = 3 one has 

a-^{q*) = an + 0.3632a^ + l.OUal + (28) 

where the use of the scale q* = 3.41/a, which is a 
sort of average value of g. In q*^ = In q^ occurring 
in the perturbative evaluation, is motivated in ref. 
[56]. The results (scale errors ommitted here) ob- 
tained from the above equation at /3 = 6.5 am- 
putated at 1-loop and then at 2-loops (plus evo- 
lution to q = 37/ro) are shown as the 4'th, 5'th 
entries in Table 3 respectively. 

^This is because the dominant contributions to the 1-loop 
coefficient df^ (s) come from tadpole diagrams and that 
these cancel when ao is replaced by ao 



The NRQCD group [29] introduce ap through 

«□= ap - (1.185 0.070n/)a|, (29) 

where the coefficient multiplying n j also depends 
on the fermion action used (Wilson in this case), 
ap is not a physical coupling; but using mea- 
surements of P and a physical scale at the same 
bare coupling one operationally assigns a value to 
a-p[q*) via eq.(29). 

The coefficients in eq.(29) are chosen s.t. 

av(g*) = ap -FO(a|). (30) 

To get to the MS-scheme one then often invokes 

IN 

"Ms(e"'^'9) = "v(g) + ^av(g)' + • • • , (31) 

where the scale is chosen so that the equation 
remains true for arbitrary nj [57]. 

Let us see how this prescription works out in 
the pure gauge theory. For example for SU(3) 
at /3 = 6.0 the scale ro/a = 5.44(26) [17]. Now 
use eq.(29) and the measurements of «□ to get an 
estimate for a\[q*) and then use eq.(31) to get 
aj^[e~'' I ^' q*) . Finally evolve via the RG with 
perturbative 3-loop /3-function to the desired en- 
ergy q. For q = 37/ro the result of this pro- 
cedure is the 5'th entry in Table 3. Note only 
one error is shown corresponding mainly to the 
scale error. The second is an unknown system- 
atic error because, one is working here at one 
value of the bare coupling, and hence an extrapo- 
lation to the continuum limit has not been made. 
The best we can do is to repeat the procedure 
for larger /3 values and observe how the results 
converge. For example for /3 = 6.4,6.5 one has 
r-o/a = 9.90(54), 11.23(21) from [17], [15] resply., 
they give the same final result, the last entry in 
Table 3. Despite the fact that with this proce- 
dure a reliable error estimate is difficult to make, 
the results seem indeed in good agreement with 
those of SF. 

7. EXTENSIONS TO FERMIONS 

So far there has been no lattice investigation 
of the actual running of a coupling in QCD. But 
encouraged by the success of the FNAL method 
to get (Xj^ in pure gauge theories (see the previ- 
ous section) , many groups have been applying the 
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method to fermions using quarkonium mass split- 
tings to set the scale. Measurements have been 
made in the quenched approximation and for dy- 
namical quarks with nj = 2 with both Wilson 
and staggered fermion actions. To obtain values 
for ftf^ in the realistic case there remain the prin- 
ciple problems of 1) the extrapolation io nj =3 
and 2) the evolution over the heavy quark thresh- 
olds (a discussion of this is postponed to sect. 9). 

7.1. Quenched approximation 

Many groups have investigated the quenched 
approximation: FNAL [21,22], NRQCD [28,29], 
UK(NR)QCD [32], KEK [31], CDHW [30] and 
ADHLM [33]. The groups all use some non- 
relativistic approximation, an expansion in pow- 
ers of V (or 1/m), where v is the velocity of the 
heavy quark in the quarkonium. However there 
are some differences in the the details of the ap- 
proaches e.g. in the number of terms which are 
kept in the expansion, the way the light quarks 
are treated, and in some cases improved actions 
are used (see e.g. [33]) for the gauge theories and 
for the fermions. Despite these differences, which 
imply different cutoff effects, there is in general 
good agreement between the results of the groups 
for both charm and upsilon (see e.g. Figs. 2 and 
3 in ref.[25]). For a given onium the relative low- 
level splittings are in good agreement with exper- 
iment. For higher levels the agreement becomes a 
bit worse e.g. there is a 4-sigma difference in the 
measured value {IP - 15')/(25' - 15) = 0.71 in 
upsilon compared to the experimental value 0.78. 

Compared to the situation reported by Michael 
last year [16], I am only aware of further progress 
by the FNAL [25] and NRQCD [27] groups, who 
have increased statistics and have concentrated 
on studying some of the systematic error sources. 
For example the FNAL group [24] found small 
cutoff effects in the IP — IS splittings, which 
translates to only a very slight a-dependence in 
Q;i^(5GeV) determined from hh IP-IS splitting. 

Setting the scale, in the quenched approxima- 
tion, using quantities with different characteristic 
scales q* , one ends up with different estimates for 
the coupling e.g. there is a 5-6 sigma discrepancy 
in the NRQCD evaluation [27] of ap(8.2GeV) be- 
tween upsilon and charm (see the first entries in 



Table 4 

Values of ap(8.2GeV) from NRQCD, nj = 3* 
corresponds to extrapolation of to Uf = 3 





scale 




4"^^(8.2GeV) 





cha. 


IS-IP 


0.1480(13) 





ups. 


IS-IP 


0.1551(11) 





ups. 


1S-2S 


0.1505(09) 


2 


cha. 


IS-IP 


0.1758(36) 


2 


ups. 


IS-IP 


0.1793(11) 


2 


ups. 


1S-2S 


0.1777(23) 


3* 


cha. 


IS-IP 


0.1940(67) 


3* 


ups. 


IS-IP 


0.1945(30) 


3* 


ups. 


1S-2S 


0.1954(43) 



Table 4). This is however not surprising because 
the quenched approximation should not describe 
the running between the different q* 's correctly, 
since is not the physical value. 

A heuristic way which was used [21] to obtain 
an estimate for nf = 4 from a quenched measure- 
ment, say in charmonium at q* ~ 5GeV, is to 
perturbatively run with nf = down to a typ- 
ical charmonium scale ~ IGeV, and then inte- 
grate back up with nf = 4. Within the quenched 
approximation alone it is however difficult to es- 
timate the systematic errors involved in such a 
procedure. 

In a recent review talk El-Khadra ^ summa- 
rized the following estimates for various errors 
contributing to the determination of a:^(6 — 
8GeV) in quenched calculations: statistical er- 
ror 1 — 3%, finite a I — 3%, perturbative 5%, sea 
quarks 5 — 8%. This gives a total of ~ 8 — 10% 
which at niz amounts to 5 — 7%. 

7.2. Dynamical quarks 

The KEK group [31] were the first to report 
dynamical quark simulations measuring IP- IS 
splittings and the plaquette expectation value P. 
They used nj = 2 Wilson quarks and extrapo- 
lated P to niq = 0. They observed that inte- 
grating back the nj = 0, 2 points, using the RG, 
the curves indeed come close at £" ~ 0.5GeV, 
which supports the original procedure adopted in 

^in Pittsburgh, June 1995 
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ref.[21]. Encouraged by this they extended the 
method to go to rij^ =3 from their measurement 
at nf = 2. Then evolving over the heavy quark 
thresholds and running to the Z-peak they ob- 
tained the result a^(mz) = 0.111(5). 

NRQCD has increased statistics on dynamical 
simulations compared to a year ago. A detailed 
report of their results can be found in these pro- 
ceedings [27]. In particular they have studied the 

sensitivity of ap"^'' to: 1) the dynamical light 
quark mass niq and 2) on the physical quantity 
used to determine (e.g. charmonium versus 
upsilon splittings). Concerning point 1) their new 
data at aniq = 0.025 is consistent with the state- 
ment that extrapolation from the bare mass value 
arUq = 0.01 used previously {a~^ = 2.4GeV) to 
realistic light quark masses has negligible effect 
compared to statistical errors. 

Concerning point 2) they find that estimates 
of the coupling ap(8.2GeV) still depend on the 
physical quantity used to set the scale, but 
the differences are not as pronounced as in the 
quenched case. The estimates from various phys- 
ical quantities have different dependences on nf 
and a comparison of nf =0 and nf = 2 results 
suggests that ap will become insensitive to which 
physical quantity is used to set the scale, once 
the correct number of dynamical fiavors nf = 3 
is incorporated. Linearly extrapolating in nf 
their values of Up^ at nf = 0,2tonj^ = 3, 
NRQCD [27] obtain the last three entries in Ta- 
ble 4 (see also Fig. 5 of their proceedings contri- 
bution [27]), which are certainly consistent with 
one another. Such observations also explain most 
of last year's discrepancy between the central val- 
ues of = 0.108(6) by Wingate et al [30] and 
0.115(2) by NRQCD [28] who used charmonium 
respectively bottomium splittings to set the scale. 

The measured value of ap has not changed dur- 
ing the past year, and NRQCD [27] still quote a 
final result of aj^{mz) = 0.115(2), where the 
main error is considered as originating from the 
present ignorance of the 2-loop coefficient for 
nf ^ 0. If this happened to be close to its rij^ = 
value (of 0.96) then their central value would in- 
crease to 0.117. 



8. THE a-COLLABORATION PROJECT 

The Alpha-collaboration has now completed 
detailed studies in pure gauge theories. Our pro- 
gram has so far been rather limited in physical 
extent; the stress has been more on precision (the 
control of systematic errors) and on the develop- 
ment of reliable methods applicable in more gen- 
eral situations. 

We are now turning to QCD. The section at 
Rome II have been investigating "bermions" [58], 
a method first suggested in ref.[59] to take fermion 
effects approximately into account by extrapola- 
tion from negative n f . 

The main next goal is to measure the SF cou- 
pling in full QCD. The analytic framework is on 
a firm footing thanks to the work of Sint [60] who 
discussed the appropriate boundary conditions 
for the fermions and studied the renormalizabil- 
ity. Sint and Sommer [61] have completed various 
1-loop computations of the SF-coupling including 
fermions; for N = 3 they obtained the term cc nf 
in Table 1. They computed with Wilson fermions 
with and without the Sheikholeslami-Wohlert [62] 
improvement term, and showed that 0(a) cutoff 
terms at 1-loop could indeed be removed (with 
SW) by a modification of the weights of plaque- 
ttes touching the boundary. As a by-product they 
also obtained the ratio of the lattice A-parameters 
for the theories with and without improvement 

Asw/Aw = exp{0.02479675n//(2&o)}. (32) 

The program development for simulations on 
the QUADRICS (APE) is well under way for full 
QCD (N = 3). I would like to draw your at- 
tention to two useful studies which have arisen 
from the preparation of this project. The first by 
Jansen and Liu [63] who have made a study of the 
Kramer's algorithm first suggested by Horowitz 
[64]. The second is an efficient method for the 
computation of the eigenvalues of the lattice 
Dirac operator, by Kalkreuter and Simma [65]. 
Tests still remain to be made as to which fermion 
algorithm should be used. Recently there have 
been some improvements [66,67] to Liischer's al- 
gorithm [68], so this looks rather promising. 

The expansion of asF along the lines explained 
in ref.[37] and the extension of the background 
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field method for computing ® in powers of 

ao to 2-loop order in QCD with Wilson quarks is 
straightforward. Such computations would cer- 
tainly be welcome, but it may take a while until 
the precision on asF which can be reached in nu- 
merical simulations of full QCD is comparable to 
the magnitude of a typical 2-loop correction. 



9. HEAVY QUARK THRESHOLDS 

Among the first to consider how to take the 
heavy quark thresholds into account were Georgi 
and Politzer [69] already in 1975. Their consid- 
erations rested on the pioneering paper of Ap- 
pelquist and Carazzone [70] in which their "de- 
coupling theorem" was stated: In a renormaUz- 
ahle field theory the heavy fields decouple at low 
momenta except for their contribution to renor- 
malization effects. 

In the MS-scheme the /3-function governing the 
running of the coupling gMS is independent of the 
quark mass. This is not in confiict with the the- 
orem since when a Green function is expressed 
as a series in the MS-coupling, the terms of the 
series are such that they convert into a series in 
an effective coupling which behaves as if only the 
light degrees of freedom are present. The explicit 
demonstration of this to 2-loops is contained in a 
series of papers by Ovrut and Schnitzer [71] and 
by Wetzel [72]. Subsequently Bernreuther and 
Wetzel [73] , assuming the AC theorem to hold for 
a momentum subtraction scheme, discussed the 
general relation (to all orders PT) between the 
couplings of the full theory and the effective light 
theory (valid in the momentum region ^ Wheavy) 
in the MS-scheme. 

How to take the heavy quark thresholds into 
account is thus quite well understood; a recent 
account has been given by Rodrigo and Santa- 
maria [74]. The discussion in principle replaces 
ad hoc phenomenological matching relations such 

as ^/Ms('^"^f') ~ ^'Ms('^"^f') (f'-"' some c ~ 2), but 
in practice such matching conditions yield results 
very close to those of a full analysis. 



^This computation is being contemplated by R. Frezzotti 



9.1. Mass dependence of the SF-coupling 

The /3-functions of all the couplings specified in 
sect. 2, are quark mass m dependent, and the run- 
ning shows the decoupling of the heavy fermions. 
This can already be seen at 1-loop; the 1-loop 
coefficient takes the form, 

ho{nf , z) = hofi + nfho,i{z) (33) 

where z = m/q (q the momentum with which the 
coupling runs) . The mass dependence of ho{nf , z) 
in the MOM-scheme was considered in [69], and 
that of the V-scheme is practically the same [75]. 




0.01 0.1 1 10 100 

z = m L 



Figure 3. Contribution of one quark to the 1- 
loop l3- function, for 9 = (dotted line), 6 = 
7r/5 (full line) for the V-scheme (dashed line). 6 
corresponds to different fermion spatial bc's. 

In their recent paper Sint and Sommer [61] ob- 
tained &o,i(-z) in the SF-scheme; it is plotted along 
with that of the V-scheme in Fig. 3. The tran- 
sition from an effectively massless quark to an 
approximately decoupled heavy quark in the SF- 
scheme is not as rapid as in the V-scheme. This 
is because in the SF-scheme there are corrections 
odd in m (cx z as z — 0, cx 1/z as z — cxd), due 
to the 'temporal' fermion bc's. As discussed in 
ref.[61] the 1/z term in the /3-function seems to 
pose a problem for the full non-perturbative com- 
putation along the lines of ref.[10], since quarks 
much heavier than l/L induce lattice artifacts 
that dominate over their physical effect if one is 
limited to lattice sizes say L/a < 20. However 
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one can in practice omit quarks with z > Zcut ~ 2, 
since if this is done then the error induced in l/(/^ 
(at 1-loop) is only about 0.003 [61], which is a fac- 
tor 10 less than the experimental error at LEP. 

10. CONCLUSIONS 

The various lattice determinations of ay^(mz) 
in pure gauge theory are consistent. These in- 
clude treatments where extrapolation to the con- 
tinuum limit is made carefully [6] -[13] and more 
direct approaches [21] where this is not the case. 

The Alpha-collaboration have initiated a pro- 
gram to study the running coupling in full QCD. 
Until results of these (or analogous) studies are 
available, in order to make statements relevant 
to present phenomenology the "shortcut method" 
which worked well in the pure gauge case has been 
applied. This has been done for quenched [25] 
and for dynamical fermions. The mq,nf,l3 de- 
pendence has been explored to some extent by the 
NRQCD group [27], but it still needs more work. 
It is crucial that the error sources are controlled 
and in particular systematic errors not underes- 
timated. Probably a realistic goal in the near 
future is to match the experimental error < 5%. 

I thank all the members of the Alpha- 
collaboration for a very pleasant joint venture. 
I am very grateful to Aida El-Khadra for send- 
ing me useful information before my talk, and to 
Junko Shigemitsu for sending a summary of the 
recent NRQCD investigations. Finally I thank 
Chris Michael - much of this review was presented 
in his talk at Cortona in February this year. 
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